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Abstract: By using Routh approximation method and related reduction method 
to deal with the model reduction of high-order systems, it is important to pre- 
serve the stability. This note proposes an effective algorithm to modify those 
methods given before. Our method is not conservative when we reduce a stable 
high-order interval system to a stable lower-order interval system. This is illus- 
trated by several numerical examples. 
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1 Introduction 

Among many system reduction methods, Routh-Pade approximation method is remarkable for lots 
of useful features such as the computational simplicity, matching of 2r time moments and preserving 
the stability ^' ^' ^. This method is introduced into interval systems ^' ^' ^. Hwang and Yang[7] 
show that the advantage of preserving stability of Routh approximation is invalid when it is used by 
direct truncation of Routh table of a stable higher order polynomial. Therefore, some works have 
been done to solve this problem 1°' ^. It is shown that it is possible to obtain an unstable interval 
Routh approximation for a stable original interval system. Then a revised method of direct Routh 
table truncation for model reduction of interval system is proposed in [8], and it claims that the 
resulting reduced order interval polynomial is guaranteed to be stable. But an example given by 
Yang [9] shows that it is not correct. And it is pointed in [9] that the principle reason of no preserving 
stability is that the interval arithmetic operations are irreversible. As a reply, Dolgin[10] formulates 
two additional conditions to avoid the situation of stability loss. This makes their method being not 
operative, concise and computational simplicity. Reference [11] proposes an improvement to 7 — 5 
Routh approximation technique for interval systems using the Kharitonov polynomials. Its advantage 
is that it can guarantee the stability of reduced models if the original system is stable. Its disadvantage 
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is that the proposed improved technique is based on Kharitonov robust stability theory, and the Routh 
table is constructed by extracting from the over bounded Routh table using Kharitonov polynomials 
with the optimization techniques. That is to say, the result is an optimization method which revised 
the Routh array through Kharitonov polynomials. This method is not concise and clear, but Routh 
approximation is. 

It is observed if the system has eigenvalue which is close to imaginary axis, whether a stable reduced 
order interval model can be derived from an interval system by using Routh-Pade approximation 
method has great relationship to the algorithm. In this case the element in Routh array of the interval 
polynomial is very sensitive to the coefficients of polynomial and is closely dependent on algorithm. 
Based on those points, we propose a new revised algorithm, which is concise, easy to compute and 
keep better preserving stability than those given before. 



2 Main Result 

Consider transfer function, 

G() = K i a o 1 + [ a i> a t] s + • • • + K-i^tiK^ = m 

[bo,bt] + [K,K]s + --- + [KX]8 n D(s) 
Let [a, b] and [c, d] be two intervals, and define 
Addition: [a, b] + [c, d] = [a + b, c + d] 
Subtraction: [a, b] — [c, d] = [a — d,b — c] 

Multiplication: [a, b] [c, d] = [min{ac, ad, be, bd}, max{ac, ad, be, bd}] 
Division: j^j = [a, b] [\,\], provided [c, d] 

What following is the introduction of our revised Routh approximation method for interval systems. 
Let G r (s) be the reduced interval systems of degree r, i.e. 

q (s) = ifoJol + [/f, + if2,f^ 2 + ■■■ + [/r-l. /rtlK- 1 A (2) 

r[S> [d ,d+] + [d^,d+}s+[d 2 ,d+]s 2 + ---+[dr,d+}s r D r ( S ) {) 

The Routh table for the denominator D(s) of the interval systems is given below, 
[^i.^i] KvbU Ksibts] [bT A M A ] ■■ 

[ & 2,l'^l] [ fe 2,2>^ 2 ] [ 6 2,3'^ 3 ] 

i b 3,V b tl} [ fe 3,2>^ 2 ] fe 6 ^] 

Kvb+J [6- ,6+] 

[ 6 i+i,i' 6 i+i,i] \-K+i,2,bt + i, 2 } 

Table 1 

where 
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And the above array is formed by algorithm (3), i.e. 



bij+ih-1,1 - h-2,ibi-i,j+i 



J i,3 



bi- 



(3) 



where b id € [6^-, 6+-], i > 3, 1 < j < 



n—i+3 
2 • 



Then the denominator of G r (s) is given by 

D r (s) = b r . :1 s r + 6 r „i,is r ~ 1 + 6 r , 2 s r ~ 2 + V_i, 2 s r ~ 3 + 



(4) 



Remark 1 : The numerator of G r (s) is given by matching the coefficients of the power series 
expansion of the system and the reduced model [4]. 

If we directly used the interval arithmetic to derive Routh stability array, then intervals of ar- 
ray(Table 1) given by (3) is enlarged. This point can be seen from the following points: 

1) All entries in each row(i > 3) are computed by using an interval, not a point; 

2) Interval arithmetic can not ensure the consistency of each pair of rows and columns; 

3) Interval arithmetic is sensitive to the computational algorithm, such as 

bi-2,ibi-ij+i bj_2,j+i&i-i,i — 6j_2,i6j_i J+ i 



H-2,3+1 



while they are same for four given numbers. 

In this paper, we focus on keeping the consistency of elements of each pair of rows and columns 
with interval arithmetic. 

Firstly, it should be pointed out that bij is increasing on 61-1,1 and decreasing on 6j_2,i6j_ij+i, 
where bij is given by 

/..../,.. 

(5) 



bij = b 



■i-2,j+l 



bj-2,ibj-ij+i 
k-1,1 



Let Aj,j = [A^-, A+-] and Vjj = [V^-, V^-] be the left and right intervals corresponding to bij. 
For i > 3, 1 < j < "~^ +3 , our revised Routh array is computed by 



A 



h3 



A A i^2,l A i^lJ+l 
Aj_2,j+1 



A 



h3 



i-2,3 + 1 



'i-2,3+1 



i-2,j+l 





A-_i,i 




A + 




A-_i,i 


v- 


2,l^i-l,j+l 








2,l V j-lj+l 


Vti,i 



i > 3,i + 4 



i = 4 



i > 3,i / 4 



i = A 



(6) 



(7) 



and 



Aj,j — Vjj 



(i = l,2;l<j<n) 



(8) 
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A 4J = [A^,A+.] */ A 4 -. = A+ 



(9) 



Therefore, our revised Routh array algorithm is computed row by row in the following fashion, 
which will be referred to as algorithm RRA. 

Step 1, evaluate the first two rows of Routh array by (8); 

Step 2, compute other rows by (6) and (7); 

Step 3, put bij = [A--, V+-], i > 3, 1 < j < ^p; 

Step 4, if A~j < Vfj for every interval element [A~-, V*-], then end; 

Step 5, otherwise, then turn to step 4. 

3 Examples 

Example 1: Consider the family of seven order interval polynomial [10] 

D(s) = [1, 2]/ + [9, 10]s 6 + [31, 35]s 5 + [71, 72]s 4 + [111, 112]s 3 + [109, 110]s 2 + [76, 84]s + [12, 13] (10) 

Firstly, the Routh table given by Dolgin's method in [8] is 

[1,2] [31,35] [111,112] [76,84] 

[9,10] [71,72] [109,110] [12,13] 

[19.63,23.79] [93.63,94.79] [73.95,82.11] 

[29.52,31.03] [73.07,77.64] [12.00,13.00] 

[37.95,42.39] [64.63,73.50] 

[17.68,28.94] [12.00,13.00] 

[42.22,52.82] 

Table 2 

Secondly, the Routh table given by Dolgin's method in [10] is 

[1,2] [31,35] [111,112] [76,84] 

[9,10] [71,72] [109,110] [12,13] 

[19.79,23.63] [93.86,94.56] [77.73,78.33] 

[29.93,30.62] [75.16,75.56] [12.33,12.68] 

[39.97,40.38] [68.95,69.18] 

[23.19,23.42] [12.46,12.54] 

[47.48, 47.56] 

Table 3 

Thirdly, the Routh table given by algorithm RRA is 
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[1,2] [31,35] [111,112] [76,84] 

[9,10] [71,72] [109,110] [12,13] 

[15.00,27.90] [86.56,101.10] [73.11,82.80] 

[19.47,335.12] [61.33,82.19] [12,13] 

[26.69,45.17] [60.42,74.63] 

[11.43,30.38] [12,13] 

[28.94,57.19] 

Table 4 



It can be observed from Tables 2,3,4 that, 

1) all of reduced-order interval models are stable; 

2) the interval elements in Table 4 are larger than those in Table 2 and Table 3; 

The fifth-order reduced interval polynomials yielded by [8] , [10] and algorithm RRA, respectively, are 

P 5 (s) = [19.63, 23.79]s 5 + [29.52, 31.03]s 4 + [93.63, 94.79]s 3 + [73.07, 77.64]s 2 + [73.95, 82.11]s + [12, 13] 

R 5 (s) = [19.79, 23.63]s 5 +[29.93, 30.62]s 4 +[93.86, 94.56]s 3 +[75.10, 75.56]s 2 +[77.73, 78.33]s+[12.33, 12.6S 

D 5 (s) = [15.00, 27.90]s 5 + [19.47, 35.12]s 4 + [86.56, 101.10]s 3 + [61.33, 82.19]s 2 + [73.11, 82.80]s + [12, 13] 

The robust step responses of uncertain systems , p^pj > ife(a) ' an< ^ Ds(s) are snown ™ * ne following 
figure: 
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The following two examples shown in [7] and [10] are reworked here. 
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Example 2: Consider a six order interval polynomial 
D(s) = [2, 2.5]s 6 + [76, 76.5]s 5 + [119, 119.5]s 4 + [111, 111.5]s 3 + [71, 71.5]s 2 + [31, 31.5]s + [9, 9.5] 



Routh tables derived by method of reference [10] and by algorithm RRA, respectively, are 



Routh Table of using method in [10] 


Routh Table of using algorithm RRA 


[2,2.5] [119,119.5] [71,71.5] [9,9.5] 

[76,76.5] [111,111.5] [31,31.5] 

[115.72,116.21] [70.18,70.47] [9.1,9.4] 

[64.86,65.16] [25.12,25.22] 

[25.38,25.48] [9.24,9.27] 

[1.51,1.54] 

[9.24,9.27] 


[2,2.5] [119,119.5] [71,71.5] [9,9.5] 

[76,76.5] [111,111.5] [31,31.5] 

[115.33,116.60] [69.96,70.69] [9.00,9.50] 

[64.95,65.07] [25.10,25.24] 

[25.03,25.83] [9.00,9.50] 

[0.55,2.47] 

[9.00, 9.50] 



Table 5 

Truncating one row in Table 5 results in the following two sixth-order polynomials: 
P 5 ( s ) = [76, 76.5]s 5 + [115.72, 116.21]s 4 + [111, 111.5]s 3 + [70.18, 70.47]s 2 + [31, 31.5]s + [9.1, 9.4] 
D 5 (s) = [76, 76.5]s 5 + [115.33, 116.60]s 4 + [111, 111.5]s 3 + [69.96, 70.69]s 2 + [31, 31.5]s + [9, 9.5] 




original system 

5nd order model using YD Mehtod 

5nd order model using RRA 
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Example 3: Consider a six order interval polynomial 



D(s) = [2.1, 2.6]s 6 + [76.1, 76.7]s 5 + [119.1, 119.6]s 4 + [111, 111.6]s 3 + [71.8, 72.3]s 2 + [31, 31.7]s + [9, 9.9] 
Routh tables of D(s) using DZ method in [8] and algorithm RRA, respectively, are 



Routh Table of using DZ method in [8] 


Routh Table of using algorithm RRA 


[2.1,2.6] [119.1,119.6] [71.8,72.3] [9,9.9] 

[76.1,76.6] [111,111.6] [31,31.7] 

[115.67,116.19] [70.82,71.35] [9.00,9.90] 

[63.98, 64.92] [24.48, 25.77] 

[24.48,27.32] [9.00,9.90] 

[-0.16,3.37] 

[9.00,9.90] 


[2.1,2.6] [119.1,119.6] [71.8,72.3] [9,9.9] 

[76.1,76.6] [111,111.6] [31,31.7] 

[115.29,116.56] [70.17,71.45] [9.00,9.90] 

[64.31,64.59] [25.07,25.17] 

[25.60,26.19] [9.00,9.90] 

[0.19,2.98] 

[9.00,9.90] 



Table 6 

It can be observed easily from Table 6 that the reduced-order interval polynomial of D(s) derived 
by DZ method is not stable, and derived by our algorithm RRA is stable while the original interval 
polynomial is stable. 

Remark 2 : It is observed from examples 1 and 2 that almost every interval element [b~j,b^~j] in 
the Routh table of using DZ method in [8] is smaller than that in the Routh table of using algorithm 
RRA, and observed from example 3 that the interval entry frg^ = [—0.16,3.17] in the Routh table 
of using DZ method in [8] contains zero. The two computed interval Routh arrays shown in table 6 
indicate that the method given in [8] cannot guarantee the success in generating a full interval Routh 
array, and our algorithm RRA can do. 

Remark 3 : If one polynomial of interval polynomials has complex roots near the imaginary axis, 
then the element in Routh array of the interval polynomial is very sensitive to the coefficients of 
polynomial and is closely dependent on algorithm. 

4 Conclusions 

The advantage of Routh approximation is to yield stable reduced-order model for stable original 
high-order systems. However, the advantage of preservation of the Routh approximation method is 
invalid when it is extended to the model reduction of interval systems. To overcome this, a revised 
Routh algorithm is proposed in this note for the model reduction of high-order interval systems. With 
our revised algorithm, a stable reduced-order interval model can be obtained as possible as. At the 
same time, almost every element in the Routh table is as large as possible. 
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